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Abstract. Using orthogonal projections, we investigate distance of a 
given elasticity tensor to classes of elasticity tensors exhibiting particu- 
lar material symmetries. These projections depend on the orientation of 
the elasticity tensor, hence the distance is obtained as the minimization 
of corresponding expressions with respect to the action of the orthogonal 
group. These expressions are stated in terms of the eigenvalues of both 
the given tensor and the projected one. The process of minimization is 
facilitated by the fact that, as we prove, the traces of the correspond- 
ing Voigt and dilatation tensors are invariant under these orthogonal 
projections. For isotropy, cubic symmetry and transverse isotropy, we 
formulate algorithms to find both the orientation and the eigenvalues of 
the elasticity tensor that is endowed with a particular symmetry and is 
closest to the given elasticity tensor. 



1. Introduction 



Computing distance between elasticity tensors is particularly important in 
inverse problems. For instance, we wish to know whether or not a material 
can be considered as isotropic with its anisotropic behavior accounted for 
by the errors of inversion, which stem from both measurement errors and 
the fact that the elasticity tensor itself is an idealization in the realm of 
continuum mechanics. 

A great portion of seismology, which is the motivation for our study, is per- 
formed under the assumption of isotropy due to its sufficient accuracy and 
mathematical convenience. However, to increase this accuracy, several seis- 
mologists formulated their equations without assuming isotropy; notably, 
Rudzki [20], [21] and, almost half a century later, Helbig [H]. Three decades 
after Helbig, Thomsen [26] popularized the concept of anisotropy among 
exploration seismologists by formulating a linear approximation whose ac- 
curacy stems from the fact that most natural phenomena investigated in 
seismology are only weakly anisotropic. 
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The purpose of this paper is to investigate the concept of distance in the 
space of elasticity tensors between a given tensor and classes of tensors ex- 
hibiting particular material symmetries, and to obtain the effective elasticity 
parameters and effective orientations for such symmetries. To facilitate this 
nonlinear process, we propose the use of invariants of the orthogonal projec- 
tions onto various symmetry classes of elasticity tensors. 

The concept of distance of a given elasticity tensor to a material-symmetry 
class was introduced by Gazis et al. [H] using orthogonal projections of 
elasticity tensors on the space of a particular class, and by Fedorov [TT] . 
who minimized the mean-square difference of the slowness surface. Three 
decades later, Norris [19] proved that the two approaches are equivalent to 
one another. Moakher and Norris [18] give a detailed presentation of distance 
of an elasticity tensor to all material-symmetry classes. In this work, we 
use the Euclidean distance, which is induced by the Frobenius norm. For 
different distances with a detailed discussion of their properties we refer the 
reader to Moakher [16] and Moakher and Norris |18| . 

In all approaches, the problem raised by Gazis et al. [13] remains open: the 
orthogonal projection, c sym , of a given elasticity tensor, c, onto the set of 
elasticity tensors belonging to a particular material symmetry class, Sym, 
depends on the original orientation of tensor c. It follows that the distance, 
d (c, c sym ), between an elasticity tensor and its orthogonal projection depends 
on the a priori orientation of tensor c. Therefore, the distance, d(c,Sym), 
between an elasticity tensor, c, and some symmetry class, Sym, can be 
obtained by finding the orientation of elasticity tensor c that minimizes the 
distance d(c,c sym ). There have been several attempts to solve this problem. 
Gazis et al. [13] give necessary conditions for the cubic case, such that the 
axes of a coordinate system coincide with the axes of its closest cubic tensor. 
Arts et al. [T] propose a method to find particular reference axes in which 
most offdiagonal terms are minimized. Dellinger [9] proposes a numerical 
method to find the rotation-symmetry axis that minimizes the expression for 
distance to transverse isotropy. Francois et al. [12] and Norris [19] propose 
numerical methods to find the closest projection. 

We begin this paper by invoking the elasticity tensor together with the corre- 
sponding Voigt and dilatation tensors. We view elasticity tensors as fourth- 
rank tensors in M 3 and consequently look for invariant properties under the 
action of orthogonal group 0(3). This point of view was introduced and 
investigated independently by Walpole [28], Rychlewski [22l|23], and Cowin 
and Mehrabadi [3 [8]. 

For a fixed orthonormal basis in M 3 , any other orthonormal basis can be 
given in terms of an orthogonal transformation A £ 0(3). The orthogonal 
projection, c sym (A), of an elasticity tensor, c, to a symmetry class, Sym, 
depends on the orientation, A, that is used to express the elasticity tensor. 
These orthogonal projections were considered first by Gazis et al. [13], and 
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explicitly obtained by Moakher and Norris [18] for all symmetry classes. 
The corresponding Euclidian distance, d (c, c sym (A)), between an elasticity 
tensor, c, and its orthogonal projection, c sym (A), may not represent the 
distance to a symmetry class, d(c,Syrn), since d (c, c sym (A)) > d(c,Sym). 
Therefore, we obtain the distance to a given symmetry class by minimizing 
d (c, c sym (A)), for A £ 0(3). To facilitate this minimization, we express 
d (c, c sym (A)) as function / (X^ m (A)), where X^ m (A) are the eigenvalues of 
the orthogonal projection c sym (A). 

In this work we prove that the traces of the associated Voigt and dilatation 
tensors, which are orthogonal invariants, also remain invariant under these 
orthogonal projections, and therefore certain relations between eigenvalues 
are preserved. These invariant relations are useful in simplifying function 
/ (\ S a m (A)). Since the minimizing process is taken over the compact group 
0(3), the minimum is reached for an orthogonal transformation, A e , which 
will be called the effective transformation or orientation. This allows us to 
obtain the distance of a given elasticity tensor to a particular symmetry class, 
and to find the orientation and eigenvalues - the effective elasticity param- 
eters - of the closest elasticity tensor that belongs to that class. We discuss 
explicitly distance to isotropy, cubic symmetry and transverse isotropy, and 
provide numerical examples for our algorithms. 

For the case of isotropy, we show that the orthogonal projection proposed by 
Moakher and Norris [18] can also be obtained using the harmonic decompo- 
sition of an elasticity tensor, studied by Backus [2] and Cowin [6]. In the case 
of isotropy, the orthogonal projection, c i— ► c* so , is invariant under the action 
of the orthogonal group. Consequently, the eigenvalues of the isotropic pro- 
jection are orthogonal invariants and therefore they can be referred to as the 
effective isotropic invariants. It implies that / (X\ so (A) , X™ (A)} is constant. 
Therefore any orientation is an effective orientation, d (c, c lso ) = d(c,Iso), 
and no minimization algorithm is required. We express this distance in terms 
of the effective isotropic parameters, as well as in terms of traces of Voigt 
and dilatation tensors for the isotropic projection. 

For the case of cubic symmetry, we use the orthogonal projection proposed 
by Gazis et al. [H] and by Moakher and Norris [H]. For an arbitrary orien- 
tation, A, of an elasticity tensor, c, its orthogonal projection, c cube (A), has 
three distinct eigenvalues, Af be {A), X c ^ be {A), Xf be {A), with multiplicities 1, 
2 and 3, respectively. The invariance of the traces of the Voigt and dilatation 
tensors under orthogonal projection, c h- > c cube , implies that the first eigen- 
value of the cubic projection is an orthogonal invariant and the other two 
eigenvalues satisfy a linear relation, which is also an orthogonal invariant. 
Therefore, the distance, d (c, c cube (A)) , can be expressed as a function of one 
eigenvalue only, / (A™ 66 ^)) , which the minimizing algorithm has to take 
into account to obtain distance d (c, Cube) . 
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For the case of transverse isotropy, we use the orthogonal projection proposed 
by Moakher and Norris |18j . For an arbitrary orientation A of an elasticity 
tensor, c, its orthogonal projection, c TI (A), has four distinct eigenvalues, 
Xj^A), X^ T (A), Xj^A), Xj^A), with multiplicities 1, 1, 2 and 2, respec- 
tively. The invariance of the trace of Voigt tensor under the orthogonal pro- 
jection c i — > c TI implies a linear relation among the four eigenvalues of the 
transversely isotropic projection. Therefore, the distance, d (c, c TI (A)^ , can 
be expressed as a function of three eigenvalues only, / [Xf 1 {A), X^ 1 (A), X^ 1 (A)^j , 
which the minimizing algorithm has to take into account to obtain distance 
d(c,TI). 

At the end of the paper, we illustrate the presented approach using examples 
of strongly and weakly anisotropic solids. For these examples, the algorithm 
provides the distances to isotropic, cubic and transverse isotropic symmetry, 
as well as the effective orientation and effective elasticity parameters for each 
of the considered symmetry classes. 

2. Elasticity tensor 

In this section, we present the notions that allow us to formulate and discuss 
distance between a given elasticity tensor and material-symmetry classes. 

The elasticity tensor is a fourth-rank tensor in M 3 , which according to Hooke's 
law relates the stress and strain tensors, 

(2.1) (Jij = CijkiSki- 

In this paper, Roman indices, i,j,..., run from one to three, and Greek 
indices, a, (3, . . ., from one to six; summation over repeated indices is implied. 
Both stress and strain are symmetric second-rank tensors in M 3 . We denote 
the space of second-rank symmetric tensors in M 3 by L 2)S (M 3 ) ; this is a 
six-dimensional linear space in which we consider the scalar product given 
by 

(2.2) e x • e 2 := Tr (ex e* 2 ) , 

where Tr stands for trace and t for transpose. Therefore, one can view 
an elasticity tensor as a linear map c : L2,s (l^ 3 ) — ► ^2,s (^ 3 ) that is 
symmetric, 

cei ■ e 2 = ex ■ ce 2 , Vei, e 2 G L 2 , s (M 3 ) , 
and positive-definite, 

ce ■ e > 0, Ve e L 2 , s (M 3 ) , e ^ 0. 

Since c is a fourth-rank tensor in M 3 , its components with respect to or- 
thonormal basis {ex, e 2 , e^} are 

Cijki = c{ei® ej) • (e fc ® ej) , 
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which satisfy the intrinsic symmetries given by Cijki = Cjiki = CkUj- For 
subsequent use, we invoke the Voigt and dilatation tensors, which are second- 
rank symmetric tensors in M 3 , given by 

(2.3) Vij = Cikjk and D^j = Cij^k, 

and discussed by Cowin and Mehrabadi [7] and Helbig |15j . 

Let us consider the group of orthogonal transformations, O (3), in the three- 
dimensional Euclidean space, R 3 . Each such transformation is linear and 
preserves the scalar product. For a fixed orthonormal basis B = {e\, e2, 63} 
in M 3 , any other orthonormal basis B' = {/1, /2, fs] can be identified with an 
orthogonal transformation A G 0(3), such that Aei = /j, for all i G {1,2,3}. 
Also, we will refer to orthogonal transformation A as an orthonormal basis. 

An orthogonal transformation A G O (3) acts on a second-rank symmetric 
tensor, e G L2 )S (M 3 ) , as 

(A o e) (n, v) := e (Au, At?) , Vu, i> G M 3 . 

An orthogonal transformation i £ O (3) acts on the elasticity tensor as 

(4 * c) (e) = A o (c (A o e )) , Ve G L 2 , s (E 3 ) . 

For the orthonormal basis B = {ei, e2, 63}, the components of ^4 * c are 

We remark here that these are also the components of elasticity tensor c with 
respect to the orthonormal basis A. 

The orthonormal basis {ei,e 2 , 63} of M 3 induces an orthonormal basis of 
L 2 , s (M 3 ) , given by 

1 

(2.4) s a (ij) = 2 2 -*<i (a <g> ej + <g) e*) , 
where 

(2.5) a : {(t, j) , 1 < i < j < 3} — > {1, 2, . . . , 6} , 

with a (i, j) = i5ij + (1 — 6ij) (9 — i — j) and Sij being the Kronecker delta. 
The components £jj = e(ej,ej) of a second-rank symmetric tensor with 
respect to the orthonormal basis {ei,e 2 ,e3} of M 3 can be arranged into a 
six-dimensional vector, 

I £22, £33) V2e 23 , V2e 13 , V2e u ) , 

with respect to basis (|2.4h . Since Lord Kelvin [24l [25], p. 110] suggested such 
an approach — albeit without using the tensorial notation, which was not 
known at the time - we refer to entities expressed with respect to basis (12. 4p 
as Kelvin's notation. An elasticity tensor, which is a second-rank positive- 
definite symmetric tensor in L2,s (R 3 ) , has the components with respect to 
basis (|2.4I) given by 

(2.6) C a p = Ce a • Ep. 
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C a/ 3 are the entries of a 6 x 6 matrix, which we write explicitly as 





Cmi 


C1122 


C1133 


\Z2CH23 


\/2cni3 


\/2c m2 " 




C1122 


C2222 


C2233 


\Z2C2223 


\Z2C2213 


\Z2C2212 


(2.7) 


C1133 


C2233 


C3333 


\Z2C3323 


\/2C33i3 


\/2C33i2 


V2CH23 


\Z2C2223 


\Z2C3323 


2C2323 


2C2313 


2C2312 




\/2cni3 


\Z2C2213 


\Z2C3313 


2C2313 


2C1313 


2C1312 




. \/2cni2 


\/2C2212 


\Z2C3312 


2C2312 


2C1312 


2C1212 . 



This formulation was introduced and investigated independently by Walpole 
[28], Rychlewski [22j [23], and Cowin and Mehrabadi [3 [8]. We refer to 
matrix (12,71) . as Kelvin's notation [2S1 P- 110]. This notation has been used 
by several researchers; notably, Fedorov [11], Helbig [15], Chapman [5] and 
Bona et al. [4]. 

For subsequent use, we invoke the scalar product given by 
(2-8) c • c := c i jkic' i j kl = C a) 3C' al3 , 

which results in the Frobenius norm given by 

(2-9) l|c|| 2 = E4^ = E C '/3- 

ijkl a/3 

We remark here that Kelvin's notation is essential for the second equality in 
both expressions (|2.8p and (|2.9p . 

The elasticity tensor's eigenvalues are positive real numbers due to its sym- 
metry and positive definiteness. We denote these eigenvalues by Ai, A2, • • • , A r 
and their multiplicities by mi,m2, . . . , m r , where m\ + 1JI2 + • • • + m r = 6. 
The norm of the elasticity tensor can be expressed in terms of eigenvalues as 

(2.10) ||c|| 2 = miA? + m 2 \l H Vm r \ 2 r . 

Scalar product (|2.8p and Frobenius norm (12. 9p are invariant under the action 
of orthogonal group O (3) . 

3. Distance to isotropy 

The harmonic decomposition of an elasticity tensor, which results in its 
isotropic and anisotropic parts, has been used by several researchers, no- 
tably Backus [2], Cowin [6] and Baerheim [3], to study properties of this 
tensor. Using this decomposition, we show that, for an arbitrary elasticity 
tensor, the isotropic part is its closest isotropic tensor, in the Euclidean sense. 
The norm of the anisotropic part, which is the component orthogonal to the 
isotropic one, is the distance between an arbitrary elasticity tensor and the 
class of isotropic tensors, as discussed by Gazis et al. |13| . Fedorov p~l] and, 
recently, by Moakher and Norris [18j. Combining these results, we discuss 
the distance to isotropy using traces of the Voigt and dilatation tensors, 
the Lame parameters, as well as the two distinct eigenvalues of the closest 
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isotropic tensor. We show that these eigenvalues are the effective isotropic 
parameters introduced by Voigt [27] and discussed by Cowin [6]. Therefore, 
the distance of an elasticity tensor to isotropy measures the deviation of its 
eigenvalues from its effective isotropic parameters. 

To derive and discuss the expression for distance, we consider an elastic- 
ity tensor whose components with respect to a fixed orthonormal basis, 
{ei,e2,es}, are c%ju. In terms of these components we consider the fol- 
lowing parameters introduced by Backus [2] and Cowin [6j for the harmonic 
decomposition of the elasticity tensor: 

(3.1) A = ^ [cim +C2222 + C3333 +4(di22 + C1133 + C2233) 

-2 (C1212 + C1313 + C2323)] , 

(3.2) fi = — [cnii + C2222 + C3333 — (cn22 + C1133 + C2233) 

15 

+3 (C1212 + C1313 + C2323)] • 

The first parameter is related to the parameters obtained by Moakher and 
Norris [18] by A = k — 2/x/3; the second one is the same as the parameter 
obtained by Moakher and Norris [18J. Using definitions (|2.3p . we express A 
and /i in terms of traces of the Voigt and dilatation tensors, 



(3.3) A = — (2TrL> - TtV) and /j = — {3TrV - TrD) , 
15 30 

which shows that these parameters are invariants under the action of the 

orthogonal group 0(3). 

According to Gazis et al. [T3], Corollary 2.1], it follows that tensor 



(3.4) 



A + 2/i 
A 
A 






A 

A + 2/i 
A 






A 
A 

A + 2/J, 










2/i 










2// 










2/x 



is positive-definite, and therefore one can view A and \i as the Lame pa- 
rameters of this isotropic tensor. Tensor c lso has two distinct eigenvalues, 
namely, 

(3.5) X\ so = 3A + 2/x = -TrD 



and 
(3.6) 



2\i 



i(3W 



TrD), 



with multiplicities mi = 1 and m% = 5. These two eigenvalues coincide with 
the two effective isotropic moduli used by Voigt |27J and Cowin |6j. 
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By a straightforward calculation we can show that (c — c lso ) 

J ) j_ c iso_ j n v j ew f this orthogonality, we denote 



(3.7) cT 
which implies that 

(3.8) ||c|| 2 = ||c iso || 2 + ||ci so || 2 . 

Therefore, the square of the distance between an arbitrary elasticity tensor 
and its orthogonal projection is 

2 /„ „iso\ ||Jsoi|2 ll„l|2 ||„isoi|2 



(3.9) d 2 (c,c iso ) = ||ci s 



c — c 



From expressions (|3.5j) and (13. 6h as well as matrix (13.41) . it follows that the 
orthogonal projection c i— > c lso is invariant under the action of the orthogonal 
group 0(3). Therefore, no minimization is necessary to obtain the distance 
to isotropy, and hence 

(3.10) d 2 (c, Jao) = d 2 (c, c iso ) = ||ci s l 2 = ||c|| 2 - ||cH| 2 - 

It follows that tensor (13. 4p is the closest isotropic tensor to a given elasticity 
tensor. Distance (13 . 10|) was obtained in terms of isotropic parameters k and 
fx by Moakher and Norris [18]. Using expression (|3.4p . we can write the 
distance (|3.10l) in terms of parameters (|3.ip and (|3.2p . 

(3.11) d 2 (clso) = \\c\\ 2 - 3 [(A + 2fi) 2 + 2 (A 2 + 2^ 2 ) 

Equivalently, using expressions (13.3(1 . we can write distance (13 . 10|) in terms 
of traces of the Voigt and dilatation tensors, 

(3.12) d 2 (clso) = llcll 2 - — (2Tr 2 L> + 3Tr 2 F - 2TryTrD) . 

15 

Also, using expression (|2.10p . we can write this distance in terms of the 
eigenvalues of the elasticity tensor and the two eigenvalues of the closest 
isotropic tensor, which are expressions (|3.5p and (|3.6|) . to obtain 

6 

(3.13) d 2 (c,lso) = J> 2 -((AD 2 + 5(Ar) 2 

a=l 

6 



(A 2 -(Ar) 2 ) + E( A '-( A "°) 2 

a=2 



Examining expression (|3- 13H . we see that an elasticity tensor is close to being 
isotropic if and only if one of its eigenvalues is close to Xf and the other five 
are close to each other and to A^ . In other words, the distance to isotropy 
measures the deviation of the eigenvalues of the elasticity tensor from the 
effective isotropic parameters. 
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4. Distance to cubic symmetry 



In this section, we discuss the distance of an elasticity tensor to the class 
of cubic tensors. We fix B = {e±,e2,es}j which is an orthonormal ba- 
sis of M 3 . For an orthogonal transformation A £ 0(3), consider AB := 
{Aei, Ae2, Ae^}, which is the corresponding orthonormal basis, and can be 
identified with the orthogonal transformation A. 

We consider the components of elasticity tensor, Cijku with respect to or- 
thonormal basis A. The following cubic tensor, which has this orthonormal 
basis as its natural basis, has been proposed by Moakher and Norris \W 









(4.1) 

where 
(4.2) 



cube 



(A) 



„cube 


„cube 


„cube 








c llll 


c 1122 


c 1122 


„cube 


„cube 


„cube 








c 1122 


L llll 


c 1122 


„cube 


„cube 


„cube 








c 1122 


c 1122 


c llll 











n cube 
zc 1212 

















n cube 
ZC 1212 


















r) cube 
zc 1212 



c mi = g ( c llll + c 2222 + C3333) 



(4.3) 

and 

(4.4) 



c 1122 — g ( c 1122 + C1133 + C2233) 



c 1212 — o ( c 1212 + C1313 + C2323) 



1 



According to Gazis et al. [13, Corollary 2.1], it follows that tensor c cube (A) 
given by expression (14.11) is positive-definite and hence represents an elastic- 
ity tensor. This elasticity tensor has three distinct eigenvalues: 



(4.5) 

(4.6) 

and 

(4.7) 



ycubt 



\ cube 
A 2 



(A) = 
(A) 



cube i q cube 
c llll ' zc 1122' 



cube cube 
c llll — c 1122 



\ cube 
A 3 



(A) 



o cube 
zc 1212) 



with multiplicities mi = 1, rn-2 = 2 and 7713 = 3, as expected in view of the 
coordinate-free characterization formulated by Bona et al. [4]. As we will see 
later, eigenvalue (14. 5j) is an orthogonal invariant. The other two eigenvalues, 
(|4.6|) and (|4.7p . are not invariant under rotations of c, and hence depend on 
orientation A. These eigenvalues were considered also by Moakher [17] and 
coincide with parameters a, c and 6, respectively, obtained by Moakher and 
Norris [18]. 
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Using matrices (|2.7|l and (|4.ip as well as expressions (14.21) - (14. 4|l . one can 
show by direct calculations that (c — c cnfee (A)) • c cube (A) = 0, which justifies 
the notation c c ± be (A) := c — c cube (A), and implies that 

(4.8) ||c|| 2 = ||c cnfee (A)|| 2 + \\cf be (A)\\ 2 . 

Therefore, the square of the distance between a given elasticity tensor, c, 
and its orthogonal projection, c cube (A), is 

(4.9) d 2 (c,c cube (A)) = \\cf be (A)\\ 2 = ||c|| 2 - \\c cube (A)f. 



This expression was obtained in terms of the cubic parameters a, b and c by 
Moakher and Norris [IB]. Using expression (|2.10l) . we can write d 2 (c, c cube (A)^ 
in terms of the eigenvalues of the given elasticity tensor and eigenvalues (|4.5j) . 
(|4.6|) and (14. 7D . Hence, we write 



(4.10) 



/ 2 2 2\ 

<i 2 (c, c cu6e (A)) = ^ A 2 - ( [Xf be (A^j + 2 (A^ be (,4)) + 3 (\ c 3 ube (A^j J 
«=i ^ ' 

2 - (A^(A)) 2 )+^ (a 2 - (A c 2 ^(^) 2 )+^ (a 2 - (Af e (.4)) 2 ) . 



a=2 v ' a=4 

which expresses the deviations of the eigenvalues of the elasticity tensor from 
its orthogonal projection. 

In order to obtain distance d (c, Cube) we have to minimize formula (I4.10p 
for A £ 0(3), which is equivalent to maximizing function 



/ (Af 6e (A),Af be (^),Af 6e (A) 



(4.11) (X? be {A)y + 2 (X? be {A)y + 3 (xf> e {A)' 2 

We now prove that \^ ube [A] is an orthogonal invariant and we can find an 
orthogonal invariant relation between X < 2 lbe {A) and Xf lbe (A). Therefore, the 
process of maximizing function (|4.11l) can be simplified. 

First, we notice that the two linear invariants associated with the elastic- 
ity tensor, TrU and TrD, remain invariant under orthogonal projection 
c i ^ c cube . Let us denote the Voigt and dilatation tensors associated with 
c and its orthogonal projection, c cube (A), by V, V cube (A) and D, D cube (A), 
respectively. In view of expressions (|2.3|) and (|4.2|) - (|4.4p . it follows that 

(4.12) TrU = TrU cu6e (,4) and Tr D = Tr D cube (A), VA e 0(3). 

Examining expressions (|3.5p and l|4.5j) . or the invariance of Tr D under both 
projections c i— > c cube and c i— ► c tso , we see that 

(4.13) Xf lbe {A) = Xf° = - TrD, 

3 
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which shows that isotropy and cubic symmetry share one effective elasticity 
parameter; it is an isotropic invariant. In other words, Af ube remains un- 
changed if the elasticity tensor is subject to an orthogonal transformation 
or its components are expressed in a different orthonormal basis. The other 
two eigenvalues of the cubic tensor, \2 Ube {A) and X^ ube (A), are not isotropic 
invariants. However, examining expressions (14. 6p . (14. 7p and (13. 6h . we see 
that their combination is an isotropic invariant: 

(4.14) 2\ c ^ be {A) + Z\f be {A) = 5A| so = Tr V - - Tr D. 

3 

Digressing, we note that, physically, the split of A?, 50 , which has multiplicity 
m 2 so = 5, into X 2 ube and A™ 66 , which have multiplicities m 2 ube = 2 and 
mjj = 3, is related to the distinction between the eigenvalues that govern 
pure and simple shears, respectively, for cubic continua. 

Expression (I4.14p is a consequence of the invariance for Tr V under both 
projections c i— > c cube and c i— > c tso , which can be written as 

(4.15) TrV cube {A) = \ c ? be {A) + 2\ c 2 ube {A) + 2,\f be {A) 

= Tr V iso = X\ so + 5A| S0 = Tr V, VA G 0(3). 

In order to obtain the distance to the subspace of cubic tensors, we must 
minimize the value of expression (|4.10l) or maximize the value of expression 
(|4.1ip . using (I4.13P and (|4.14|) . Since this minimization is performed over 
the compact group 0(3), it follows that the minimum of expression (I4.10p 
can be reached for an orthogonal transformation, A e . Hence, there exists 
A e S 0(3), such that 

(4.16) d (c, Cube) = d (c, c cube {A 

Therefore, c cube (A e ) is the closest cubic tensor to elasticity tensor c. Eigen- 
values Xl ube (A e ), \2 Ube (A e ), Xf lbe (A e ) represent the effective cubic parame- 
ters, while A e is the effective cubic orientation for the elasticity tensor c. 

Since changes of orthonormal basis in M 3 affect neither the six eigenvalues 
of the elasticity tensor nor the isotropic invariant, Af^ e , using expression 
(|4.14p . we state expression (|4.10l) as 

(4.17) d 2 (c,c cube {A)) = d 2 (c,lso)-™(\r-\? be (A)] 

2 /„ t"^„\ ( \iso \cube{ 



= <?(c,i8o)-f [xr-xr e (A) 

Both expressions 10 (\ 2 so - \ c 2 ube (A)f /3 and 15 (A| so - Xf be {A)f /2 repre- 
sent the square of the distance between c cube (A) and the class of isotropic 
tensors. In other words, we can rewrite expression (|4.17l) as 

(4.18) d 2 (c, c cube (A)^j = d 2 (c,lso) - d 2 (c cube (A),lso 
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Expression (I4.18P was implemented in the same context by Moakher and 
Norris [18] . 

If in expression (|4.18|) we choose A to be the effective cubic orientation A e 
and use equation fl4. 16j) . we obtain 

(4.19) d 2 (c,Cube) = d 2 (c,lso) - d 2 (c cube (A e ),lsc?j , 

which illustrates the fact that, for a given elasticity tensor, its distance to 
isotropy is greater than its distance to any other material-symmetry class. 

If in expression (|4.10|) we choose A to be the effective cubic orientation A e 
and use equation (I4.16p . we obtain 

(4.20) 

d 2 (c,Cube) = £ A*- ((\T be {A e )) 2 + 2 (\T be (A e )) 2 + 3 (a^> e (A e )) 2 ) 

a=l ^ ' 

\\ - (Ar be (A e )) 2 )+^ (\l - (xT be {A e )f\+Y^ Ul - (\? be (A e )) 

which expresses the deviations of the elasticity tensor from cubic symmetry, 
in terms of the effective cubic parameters. In other words, an elasticity 
tensor is close to having cubic symmetry if one of its eigenvalues is close to 
Af be , two of them are close to one another and to \2 Ube (A e ), and the other 
three are close to each other and to X^ ube (A e ). 



5. Distance to transverse isotropy 



In this section, we discuss the distance of an elasticity tensor to the class 
of transversely isotropic tensors. We consider the components of elasticity 
tensor, Cijki, with respect to orthonormal basis A. The following transversely 
isotropic tensor, which has this orthonormal basis as its natural basis, has 
been proposed by Moakher and Norris |18j : 



(5.1) 

where 
(5.2) 



TI 



{A) 



TI 
1111 


r TI 
L 1122 


r TI 
c 1133 








TI 
1122 


c llll 


f .TI 
c 1133 








TI 
1133 


C TI 
c 1133 


r TI 

c 3333 

















n r TI 
zc 2323 

















zc 2323 


























r TI 

c llll 



1 



-1111 



(3ciin + 3c 2 222 + 2CH22 + kc\2Yl) , 



„TI 
-1122 



(5.3) 



r TI 
L 1122 



(Cllll + C 22 22 + 6C1 



122 



4c 



1212 J 
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(5-4) cf/33 = - (di33 + C2233) , 

(5-5) C3333 = C3333 

and 

(5-6) = - (C 23 23 + C1313) • 

According to Gazis et al. [HI Corollary 2.1], it follows that tensor c TI (A) 
given by expression (|5.ip is positive-definite and hence it represents an elas- 
ticity tensor. This elasticity tensor has four distinct eigenvalues: 
(5-7) 

„TI , r TI _i_ r TI , /( TI 7 r TI _ r TI \2~Q ( r TI \ 2 

» / a\ 1111 1122 3333 y \ 1111 1122 L 3333,< ~r C H L 1133,' 

*l {A) — - , 

(5.8) 



r TI ,„TI , r TI _J( r TI , JI _ r TI \ 2 ,o( r TI \ 2 
TI _ c llll ' L 1122 c 3333 V V 1111 1122 L 3333j l L 1133^ 



(5-9) ^ 1 {A)=c{( 11 -c{( 22 
and 

(5.10) \J I (A) = 2cli 23 , 

with multiplicities mi = 1, ni2 = 1, m 3 = 2 and 7714 = 2, as expected in 
view of Bona et al. [4]. Again, unlike eigenvalues (|3.5p and (|3.6p . these 
eigenvalues are not invariant under rotations of c, and hence they explicitly 
depend on orthonormal basis A. These eigenvalues were also considered by 
Moakher [17]. The relations between these eigenvalues and the parameters 
obtained by Moakher and Norris [18] are 



a + b±J(a-bY + c 2 mr 
(5-11) Mi(A) = ^— -1 , \f {A) = /, XV {A) = g. 



Using matrices (|2.7p and (|5.ip as well as expressions (15. 2|) - (15. 6p . one can 
show by direct calculations that (c — c TI {Aj} ■ c TI (A) = 0, which justifies 
the notation c^ 7 (A) := c — c TI (A), and implies that 

l|c|| 2 = ||c"(A)|| 2 + || C ^^)|| 2 . 

Therefore the square of the distance between a given elasticity tensor, c, and 
its orthogonal projection, c TI (A), is 

(5.12) <? (c,c TI (A)) = ||cl'(,4)|| 2 = ||c|| 2 - ||c T V)|| 2 . 

Expression (I5.12P was obtained by Moakher and Norris [18] using five param- 
eters: a, 6, c, / and g. Using expression (|2.10|) . we can write d 2 (c, c TI {Aj) 
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in terms of the eigenvalues of the given elasticity tensor and the eigenvalues 
(|5.7j) - (I5.10p of the orthogonal projection. Hence, we write 

(5.13) 

d 2 (c,c T \A)) = A*-((AfV)) 2 + (A^(A)) 2 + 2 (AP(vl)) 2 + 2 (A^V)) 2 ) 

a=l 

= (a? - (xr(A)) 2 ) + (M - (a^V)) 2 )+£ (a 2 - (aF(^)) 2 )+E (xl - (aTV)) 2 ) , 

a=3 a=5 

which expresses the deviations of the eigenvalues of the elasticity tensor from 
those of the corresponding orthogonal projection. 

Using expression (13 . 13[) . we write expression fj5 . 13|) as 

(5.14) d 2 {c,c TI {Aj) = d 2 (c,lso) + (Af ) 2 + 5 (A^) 2 
- ((AfV)) 2 + (A^)) 2 + 2 (AF04)) 2 + 2 (AlV)) 2 ) • 

In order to obtain distance d (c, TI) , we have to minimize the above formula 
(|5.14p . for A in 0(3), which is equivalent to maximizing function 

(5.15) /(Af(A) ! Af(A),Af(A),Af(A)) = 
(AfV)) 2 + {XV{A)) 2 + 2 (Af(A)) 2 + 2 {XV{A))\ 

Even though these expressions are not invariant under orthogonal transfor- 
mations, the two linear invariants associated with the elasticity tensor, Tr V 
and Tr D, remain invariant. We denote the Voigt and dilatation tensors as- 
sociated with c and its projection, c TI (A), by V, V TI (A) and D, D TI (A), 
respectively. In view of expressions (12. 3|) and (15. 2|) - (|5.6p . it follows that 

(5.16) TrV = TrV TI (A) and Tr D = Tr D TI (A), \/A G 0(3), 
which are analogous to expressions (|4.12|) . 

Examining expressions (15.7(1 - (15 . 10|) , or the invariance of Tr V under pro- 
jections c i— > c TI and c i— > c 4SO , we obtain 

(5.17) TtV TI (A) = Xj I (A) + X% I (A) + 2Xj I {A) + 2Xj I (A) 

= Tr V iso = X\ so + 5Ajf° = Tr V, VA G 0(3). 

We note that Tr D TI can be expressed in terms of two eigenvalues Xf 1 , X^ 1 
and another orthogonal invariant, which is used by Bona et al. [4j for char- 
acterizing transversely isotropic symmetry but does not appear explicitly in 
the distance function; it is referred to by Bona et al. [4] as 7. Therefore, the 
invariance of Tr D TI does not have an influence on minimizing expression 

(prop . 

To obtain the distance to the class of transversely isotropic tensors, we mini- 
mize expression (15 . 13|) under the restriction given by expression (15 . 1T|) . This 
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is equivalent to maximizing function (15 . 15|) under the same restriction. Using 
expression (I5.17p . we write expression (I5.14p as 

(5.18) d 2 (c, c TI (A)) = d 2 (c,lso) + (X\ so ) 2 + 5 (Xf ) 2 
(\r(A)) 2 + (\T I (A)) 2 + 2(\J I (A)) 2 

+\ (Al so + 5A| so - (Af (A) + \l\A) + 2\l I (A)))' 

which expresses the square of the distance to transversely isotropic projec- 
tion. 

The minimum value of (I5.18P represents the square of the distance of the 
given tensor, c, to transverse isotropy. Since the minimization is performed 
over the compact group 0(3), it follows that the minimum can be reached 
for an orthogonal transformation, A e . Hence, there exists A e G 0(3), such 
that 

(5.19) d{c,TI) = d(c,c TI (A e )) . 

Therefore, c TI (A e ) represents the closest isotropic tensor to elasticity tensor 
c. Eigenvalues Xf^Ae), A^ 7 (A e ), X^^Ae), Xj^Ae) represent the effective 
transverse isotropic parameters, while orthonormal basis A e gives the effec- 
tive orientation of the closest transversely isotropic counterpart to c. 

If in expression (|5.14p we choose A to be the effective transversely isotropic 
orientation, A e , and use equation (|5.19|) . we obtain 

(5.20) d 2 {c,TI) = d 2 (c,lso) - d 2 (c T1 (A e ),lso) . 

If in expression (15 . 13|) we choose A to be the effective transversely isotropic 
orientation, A e , and use equation (|5.19l) . we obtain 

(5.21) 

d 2 (c,TI) = £ A^-((Ar^ e )) 2 + (AlVe)) 2 + 2 (AP(^ e )) 2 + 2 (Af (A e )) : 

a=l 

= (A? - (Ar04 e )) 2 ) + (A| - (xV{A e )) 2 ) + j2 (X 2 a - {Xl\A e )) 2 )+j2 (x 2 a - (A 4 "(A e )) : 

which expresses the deviations of the elasticity tensor from transverse isotropy, 
in terms of the effective transversely isotropic parameters. In other words, 
an elasticity tensor is close to having transversely isotropic symmetry if one 
of its eigenvalues is close to A^(A e ), one is close to X^ I {A e ), two of them 
are close to one another and close to Aj 7 (^4 e ), and the other two are close 
to one another and close to X^^Ae). 
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6. Examples 

In this section, we use the expressions derived above to illustrate the finding 
of the distance between a given elasticity tensor and isotropy, cubic symmetry 
and transverse isotropy. We consider both strongly and weakly anisotropic 
materials. In both cases, we present a density-scaled elasticity tensor in the 
Kelvin's notation stated in expression (12.7(1 . 

6.1. Strong anisotropy. Let us consider a strongly anisotropic case, namely, 

" 56.60 8.98 3.45 

8.98 56.60 3.45 

, s 3.45 3.45 16.43 

^ ' 3.60 

3.60 

47.62 

which corresponds to biotite and exhibits transverse isotropy, as quoted by 
Thomsen [26]. Its anisotropy parameters defined by Thomsen [26] are 7 = 
6.12, 5 = —0.39, e = 1.22; their magnitude indicates strong anisotropy. The 
eigenvalues of this tensor are 66.07 and 15.95 with multiplicity one, and 47.62 
and 3.60 with multiplicity two. Expression (16. ip is considered with respect 
to a natural basis of the transversely isotropic medium. 

6.1.1. Relation to isotropy. Using equations (|3.5p and (13. 6p . we obtain the 
eigenvalues of the closest isotropic tensor, X\ so = 53.80 and X™ = . . . = 
Xq° = 26.13, which are the effective isotropic parameters, as well as A = 9.22 
and [i = 13.06, which are the Lame parameters of the closest representative 
of the isotropic subspace. Using these values and any of expressions (13. lip 
- (I3.13p . we obtain the distance to isotropy: d(c,Iso) = 53.59 km 2 /s 2 . 

6.1.2. Relation to cubic symmetry. Using expressions (14. 5p . (14.6(1 and ((4.71) . 
we obtain the eigenvalues of the corresponding cubic tensor, with respect to 
the fixed natural basis of elasticity tensor c: Af be = 53.80, Xf> e = 37.92 
and A§" 6e = 18.27. For an arbitrary orthonormal basis A, consider c cube (A) 
the cubic orthogonal projection. Invoking one of the equalities in expression 
(|4.17p and using d(c,Iso) and A?> so from above, we write 

(6.2) d 2 (c,c cube {A)) = d 2 {caso)- l ^(xr-\T he (A)) 2 

= 53.592 _ 10 13 _ X cube^ 2 _ 

To find the distance to the cubic-symmetry subspace, we look for A™ 6e (^4) 
that minimizes this expression. The resulting distance d (c, Cube) = 49.07 km 2 /s 2 



"km 2 " 
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is achieved by rotation of the natural basis of elasticity tensor c with the ef- 
fective transformation 







0.58 







1 



-0.82 -0.58 



-0.82 



The resulting cubic orthogonal projection c cube (A e ) represents the closest 
cubic tensor to c. Its eigenvalues, which are the effective cubic parameters 
for elasticity tensor {63]), are Xl ube (A e ) = 53.80, Xf* be (A e ) = 12.64 and 
X'^ lbe (A e ) = 18.27. Minimizations described herein lend themselves to such 
software as Maple and Mathematica. 



6.1.3. Relation to transverse isotropy. Minimizing expression (15 . 18|) subject 
to Xj 1 , X^ 1 , A3 1 and Xj 1 being positive, we obtain d (c, TI) = 0, as expected. 
Since tensor (|6.1h is expressed in natural coordinates, we can obtain the 
same result using expressions (|5.7p - (I5.10P to calculate the eigenvalues of 
the effective transversely isotropic tensor, and verify that they are the same 
as the eigenvalues of tensor (16. ip . 



6.2. Weak anisotropy. Let us consider a weakly anisotropic case, namely, 



(6.3) 



4.00 


2.06 


2.10 


-0.07 


0.01 


-0.03 


2.06 


3.83 


1.96 


0.17 


-0.07 


0.18 


2.10 


1.96 


3.96 


0.16 


0.04 


-0.13 


-0.07 


0.17 


0.16 


2.00 


0.22 


-0.14 


0.01 


-0.07 


0.04 


0.22 


1.76 


0.02 


-0.03 


0.18 


-0.13 


-0.14 


0.02 


2.22 



km 
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which are the density-scaled elasticity parameters used by Dewangan and 
Grechka |10j . This is a generally anisotropic medium: it does not exhibit 
any material symmetry. The eigenvalues are distinct from each other: 8.02, 
2.39, 2.16, 1.86, 1.82, 1.52. 



6.2.1. Relation to isotropy. The eigenvalues of the closest isotropic tensor 
are Af = 8.01 and A?> so = • • • = Ag so = 1.95, which are the effective isotropic 
parameters, and the Lame parameters are A = 2.02 and fj, = 0.98. Herein, the 
anisotropy is very weak: the distance to isotropy is 0.724 km 2 /s 2 . Notably, 
by inspection of the corresponding entries in matrices (13. 4|) and (16. 3p or by 
comparison of the eigenvalues of c and c iSO , we can expect the closeness of 
this tensor to isotropy. 
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6.2.2. Relation to cubic symmetry. Proceeding in a manner analogous to the 
one described in Section [6. 1 .21 we obtain d(c,Cube) = 0.64 km 2 /s 2 , where 
the effective orientation of the closest cubic tensor to elasticity tensor (16. 3ft 
is 

0.77 0.36 0.53 



-0.60 0.13 0.79 



0.21 -0.93 0.31 j 

The effective cubic parameters of the elasticity tensor (|6.3p are Xf lbe (A e 
8.01, Xf be {A e ) = 0.59 and X^ ube (A e ) = 2.08. 



6.2.3. Relation to transverse isotropy. Proceeding in a manner analogous 
to the one described in Section [6 . 1 .31 we obtain d{c,TI) = 0.57 km 2 /s 2 , 
where the effective orientation of the closest transversely isotropic tensor to 
elasticity tensor (|6.3|) is 



A P 



0.76 0.53 0.38 



-0.58 0.29 0.76 



0.29 



-0.80 0.53 



The effective transversely isotropic parameters of the elasticity tensor (16.3 
are Xj\A e ) = 8.01, A^(^4 e ) = 1.62, X^ I (A e ) = 2.05 and Xj T (A e ) = 2.02. 
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